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The main properties of the interaction of ultra low-energy antiprotons (E < 10 6 a.u.) with 
atomic hydrogen are established. They include the elastic and inelastic cross sections and Protonium 
(Pn) formation spectrum. The inverse Auger process (Pn + e — > H + p) is taken into account 
in the framework of an unitary coupled-channels model. The annihilation cross-section is found 
to be several times smaller than the predictions made by the black sphere absorption models. A 
family of pH nearthreshold metastable states is predicited. The dependence of Protonium formation 
probability on the position of such nearthreshold S-matrix singularities is analysed. An estimation 
■ for the HH annihilation cross section is obtained. 
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& ; I. INTRODUCTION 

The unique features of the LEAR (low-energy antiproton ring) facility at CERN made recently possible the synthesis 
of few antihydrogen atoms Q]. This effort would be pursued with the antiproton deccelerator (AD) project § and 
would make possible the storage of sensible amounts of antihydrogen. This project reinforces the already active 
O • interest to investigate several theoretical and experimental problems in the physics of antimatter 

In view of storing antimatter in traps, it would be interesting to have some theoretical calculations of the rate at 
which antiprotons (p) and antihydrogen (H) annihilate with the residual gas. This process should be evaluated at 
subkelvin temperatures, the optimal energy domain for an effective synthesis and trapping of antimatter. 
. From a theoretical standpoint, a specific feature of systems like H + p or H + H, containing pairs of unlike charged 
£C) ' heavy particles (p and £>), is the possibility of rearrangement followed by Protonium (Pn) formation. Indeed even at 
zero p kinetic energy Protonium can be produced in states with principal quantum number n < 30. 

Since the pioneer work of E. Fermi and E. Teller ^| in 1947, this problem has been treated by several authors 
|p|, ^0|Jr^1 in the energy range going from a fraction to tens a.u.. The usual approach is based on the following two 
assumptions: i) the separation of electronic and nuclear motion and ii) the classical treatment of the antiproton 
— ^ j dynamics. The aim of the present work is to provide a correct description of the ultra-low energy limit, i.e. Tp < 1CP 6 
a.u., an energy domain in which the above mentioned assumptions are no longer valid. As a consequence we are 
. i-h ' definitely faced to a quantum three- or four-body problem. 

We will consider in what follows the annihilation of ultra slow antiprotons with atomic hydrogen in the framework 
pH ' of an unitary coupled channel approach. This process will be identified to the free Pn formation: 



H -> Pn* + e (1) 



as far as the direct annihilation at such low energies can be neglected and that the further evolution of Pn states will 
result into annihilation. 

In our treatment the Protonium formation as well as the virtual rearrangement process, i.e. the inverse reaction of 
(0): 

p + H -> Pn* + e^p + H (2) 

will be properly taken into account. 

We will show that the energy dependence of the inelastic reaction probability is determined by a rich spectrum 
of nearthreshold S-matrix singularities, corresponding to Hp nearthreshold metastable states generated by the long 
range charge-dipole interaction. Finally we will give an estimation of the HH annihilation cross-section in the energy 
range from 10 -8 to 10~ 4 a.u. 
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II. THE FORMALISM 



The adequate formalism for the three-body problem are the Faddeev equations |13| , according to which three 
possible asymptotic clusters have to be explicitly described. In the case of slow antiprotons scattering on hydrogen 
only two of them are physically important. The first one is the (pe)p cluster, which corresponds to the elastic 
channel H + p — > H + p. The second one is the {pp)e cluster, corresponding to the Protonium formation channels 
H +p~* Pn* + e. 

A direct solution of the Faddeev equations for such a problem is made difficult by the big number of open channels 
containing fast oscillating asymptotics. In this section we will develop a formalism which enables to take into account 
the correct behavior of the three-body wavefunction in both mentioned asymptotic clusters and benefit from the small 
value of the electron-antiproton mass ratio. We will show that this approach gives the scattering observables with an 
accuracy of the order of ~ 10% by using a very limited number of channels. In the same time it provides a transparent 
physical understanding of the low energy three-body dynamics in the pH reaction, which was the main aim of our 
study. 



A. The coordinate system 

The Jacobi coordinates for the 3-body problem are connected with the possible different asymptotic clusters. The 
coordinates corresponding to the elastic channel are defined by 



„ „ m e r e + MpRp 

XV — ixp — 



(3) 



M p 



where r e , R p and Rp are respectively the electron, proton and antiproton coordinates, and m e M p are the electron 
and proton mass. 

It turns out that this frame is also convenient for describing the Protonium production channels. Indeed, the p-p 
distance p = R p — Rp, coincides within r(m e /M p ) with R in equation (||). We will show that knowing the three-body 
wavefunction at r -C p(M p /m e ) is enough to get a good approximation of the scattering observables. Thus we can 
substitute p by R and take into account, if necessary, the difference between p and R in a perturbative way. We will 
use hereafter the unique coordinate system (^|). 



B. Three-body wavefunction 



The three-body wavefunction is represented as a sum of two components, corresponding to the two considered 
clusters: 

$(r J R)=*i(r,R) + * 2 (r,R) (4) 
Component <I>i is supposed to describe the elastic channel and can be written as 

$i(r,R)=£> a (r) %Q (R) (5) 

a 

where a is the set of quantum numbers labeling the Hydrogen atomic states <p a as well as the corresponding antiproton 
wavefunction \a- We will denote from now by channel each term in the three-body wavefunction expansion like (|J). 

For incident antiproton energy much smaller than the first Hydrogen excitation threshold, it is convenient to select 
from (|^) only the contributions which do not vanish in the asymptotics. This gives the simple form: 

$ 1 (r,R)=0 ls (r)x(R) (6) 

It is useful to introduce a projection operator P, which acts in the three-body states space and projects on the 
subspace of Hydrogen states corresponding to open channels: 

P = ^\Ki m ){<i>ni m \®i (7) 



2 



In our case the sum is limited to the Is state only. Component $1 is then written as: 



$1 = Pi> (8) 

The second component $2 describes all the remaining channels and can be written in the form: 

* 2 = (i-P)$ (9) 

This means that all the electron states except Is contribute into $2 ensuring the orthogonality of both components. 
The $2 component contains terms which correspond to the Protonium formation channels. In order to explicitly take 
into account the asymptotic behavior of the cluster (pp)e, $2 is expanded in a complete set of the (pp) eigenfunctions 
/MR): 

$ 2 (r,R)=]T 3/3 (r)MR) (10) 



where gp (r) are unknown expansion coefficients representing the electron wavefunctions in the channels characterized 
by Protonium quantum numbers f3. 

Let us remark here that at this level no approximation has been done. In particular the truncation done in the 
choice (JsJ) of $1 is balanced in the functions g@ of the second component. 

C. Equations 

The Schrodinger equation for a three-body vector state |$ > reads: 

(H ep +H^ + W^-E)\$>=0 (11) 



with: 



the Hydrogen Hamiltonian, 



the Protonium Hamiltonian and 



H ep = -J-A r -1 (12) 



1 

2M~ n |R + r(m e /M p )| 



W ex = (14) 

e P |R-r(l-m e /M p )| 1 ' 

the electron-antiproton interaction potential. M is the pH reduced mass which, neglecting m e /M p terms, will be 
hereafter approximated by the Protonium reduced mass M w M p /2, E = eb + Ep the total energy, eb the Hydrogen 
ground state energy and Ep the center of mass energy of incident antiproton. Note that all the spin degrees of freedom 
are neglected. We will first neglect the term r(m e /M p ) and substitute the exact H p p and W^p by the approximations: 

ELnri = A 7? — 

pp 2M R 

^ = WT\ 

Using (||) and (||) we obtain the following coupled equations for the components |$j >: 

{H p p + PW e pP - £ p -)|$i) + PW e p{\ - P)|$ 2 ) = 

(15) 

(H ep + Hpp + (1 - P)W e p{l - P) - £)|$ 2 ) + (1 - P)W e pP\3>i) = 
We used here the fact that P commutes with H ep and H p p. 
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The corresponding equations for the antiproton x(^) an d electron wave- functions gp{r) could be obtained by 
substituting (g) and ([l(]) into (|l^) . Due to the choice of the wavefunction components the solution of such a coupled 
equations system will correctly describe the asymptotic behavior of the three-body system. This procedure will 
however remain formal, for it leads to an infinite set of coupled channels, including the closed ones, characterized by 
a continuous pp momentum variable. To construct an equation system suitable for practical calculations, we should 
first analyze the contribution of different channels in the expansion (0). In particular, the one coming from the 
continuous spectrum. For such a purpose we represent <E>2 as a sum of two components: 

$2 = $2 + $2 (16) 

where 

$|(r, R) =F$ 2 (r, R) = £ gptffpQl) (17) 





*5(r,R) = (l-F)<h(r,R)= £ g (r)f (R) (18) 

and 



/3=n„ 



n m ax is a certain set of Coulomb quantum numbers chosen in such a way, that $2 contains all the open channels 
and, eventually, a limited number of closed ones. The sum from n max to infinity includes also the integration over 
the continuous p momentum. The $2 component describes the dynamics of Protonium formation and includes the 
corresponding asymptotics of the three-body wavefunction : 

lim $f(r,R)=V^(R)^/ l +(r) 



Here Sp are the S-matrix elements for the Protonium formation and ht(r) are the outgoing electron waves in the 
channel with quantum numbers f3. The component $2 contains only closed channels. 

At big R the component $2 vanishes due to the Coulomb bound state wave functions fp. On the contrary, the 
contribution from $2 is essential as far as it contains non vanishing terms coming from the pp states in the continuum 
corresponding to the virtual excitations and break-up. It is shown in Appendix that $2 actually describes the effect of 
the long-range polarization which can be taken into account by introducing in the elastic channel the local potential 



with a(R) ensuring that for R >> rs (H-Bohr radius) a(R) — » — ay the H dipole polarizability. The following 
approximation pi for a(R) was found to be suitable in practical calculations: 

a(R) = ^(r 5 + a d R 4 + 2a d R 3 + ^a d (R 2 + R + W e~ 2R - a d 

It is qualitatively shown in Appendix and proved by numerical calculations that at distances R « rs, $2 dominates 
and the main contribution to the wavefunction comes from the channels with n=26-40, thus we can choose V po i(R <ti 
r B ) -> 0. 

We will first consider the case of total Hp angular momentum L equal zero. The characteristic incident p energy, 
below which S-wave dominates in the elastic channel will be determined later. 

It was shown in |l5|,|l6| that in the L — case, Protonium is primarily produced in states with angular momentum 
I = and 1 = 1. The physical reason is that for the open channels with n=26-30, which dominate the reaction 
amplitude, the electron is ejected with rather small momentum fc e , and the centrifugal barrier reduces the probability 
to find a slow electron (and consequently Protonium) with high angular momentum. 
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The preceding results enable us to construct a model which, including a limited number of channels, correctly 
describes the asymptotic behavior of the three-body wavefunction. These channels dominate the reaction amplitude, 
while closed channels corresponding to the continuous spectrum of Protonium states, are taken into account by means 
of the polarization potential (20). In practical calculations we have included the Protonium channels with principal 
quantum number n=10-40 and angular momentum I = 0, 1. Numerical checks showed the stability of the results when 
increasing the number of included channels. The considered equation system has the form: 



2M dR 



L(L + 1) 



2MR 2 



V CS (R) + V pol (R) -E + s B ) X (R) 



J fa»(r)Woi(r,R)f n ,i(R)itg n ,i(r)dr = 



(21) 



E 

n'V 



-di- 



i + 



2m r 2mr 2 



+ E n 5nn',w + nU nn > t a> (r)7T 



9n',v(r) 



+ TnMr) / f nil (R)Wi (R,r)x(R) dR = 



(22) 



Here : 



V CS (R) 
W ol (r,R) 

U n n',ll' 
TT 



1+ R ]e 



-IB. 



r l + l V2TTI 



J / lM (R)/„ % r(R)^-^d 3 R 



= 1 - PS, 



0/ 



and E r , = —E 



M 
In 2 - 



D. Effective potential method 

The numerical solution of (|2l[-p2|) is a rather difficult task, as far as each of the coupled integro-differential equations 
includes fast oscillating functions \ an d fn- 

A way to overcome this difficulty is by means of the effective potential method. This approach turned to be an 
efficient tool for both a qualitative understanding and precise numerical treatment of the problem. The Hp scattering 
observables are given by the first equation for antiproton wavefunction v ( f2l| ) . The effects of the remaining coupled 
equations ( p2|) are taken into account by transforming the system ( ^l| - p2| ) into one equation for \ which contains a 
complex nonlocal effective potential: 

'Jm 9 * + L( 2Mte + Vcs(R) + fe// + Vpol{R) - E + £B ) *W = ( 23 ) 

with 

Ve ff = l/n,l><^i.l^«(r,ii)^S(ry)*W'ro(r',iJ')|0i.></nV'| (24) 



nn' ,11' 

Here f (r, r l ) is the Green-matrix of coupled equation system for electron wave-functions: 



= ( {-^di + ^ + En \ j e + *v!L>*) (25) 

The whole problem is then splitted into two parts: to calculate the effective potential and to solve the one-channel 
problem for antiproton scattering in a complex nonlocal potential. 

The benefits of such an approach are several. On one hand the Green function (25) is calculated by solving a 
coupled equation system for smooth electron wavefunctions, while the fast oscillating pp wavefunctions are explicitly 



-l 
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introduced by the well-known Coulomb states. On the other hand the effective potential (|24|) practically does not 
depend on the p incident energy in the domain Ep <C 0.01 a.u. 
(from Pn state with n=30) is about 0.02 a.u. and for p energies less than this value, the Green-matrix (|25j) is not 
sensitive to incident antiproton energy. This means that, once calculated for Ep — 0, the effective potential can be 
used in the whole energy range of interest and this radically simplifies the calculations. 

From a physical point of view it seems also more natural to analyze the properties of the Hp system in terms of a 
modified one channel problem. The main features of the effective potential as they appear from our calculations are 
the following: 

1. The imaginary part of V e ff vanishes at distance R rs 1.8re, which corresponds to the mean radius of the last 
Protonium open channel (n=30). In Fig. g Im[V e ff] (R, R') for R' = R is plotted as a function of R. 

2. The imaginary part of V e f f(R, R 1 ) is sharply peaked around its diagonal R' = R. Nevertheless for R < tb, the 
nonlocality range is of the same order as the antiproton wavefunction x(R) oscillation period. The profile of 
Im[V ef /] (R, R') for R = 0.5 is shown on Fig. |. 

3. The profile of the real part of V e ff is plotted in Fig. [| Its nonlocality range is bigger than for the imaginary 
part. It vanishes at R « 3t\b and dominates over the polarization V po i and the Coulomb screened V cs potentials 
in the range 1 < R < 3 re 



15 lq|. The minimum energy of the ejected electron 



III. RESULTS 



In this section we will present the main results obtained in the coupled-channels model and discuss the physical 
reasons of certain scattering observables behavior. 



A. Scattering observables 



The Hp complex scattering length is found to be: 

a = (-7.8 - ill.5)r B 
The corresponding elastic cross-section at zero energy is: 

a e i = 2426.4r| 

We remark the relatively big value, in the atomic scale, of the scattering length imaginary part. Such value is a 
consequence of the long range polarization forces. By switching off V po i in ( p3[ ) the value obtained is substantially 
reduced to Im(a) = 0.2tb- The capital role of the polarization forces in the low energy Hp dynamics will be discussed 
in the next subsection. 

We have calculated the energy dependence of the inelasticity S% for several partial waves. The results are shown on 
Fig. ^ for p incident energies in the range from to 10~ 6 a.u.. The inelasticity turns to be less than 0.1 for Ep < 10~ 8 
a.u. and does not become greater than 0.5 in the energy domain of interest. One can also see in this figure that the 
scattering length approximation is valid for energies less than 10~ 8 a.u. The results for I ^ have been calculated 
under the assumption that the effective potential ( p4[ ) weakly depends on total angular momentum L in the energy 
range of interest. As one can see, S-wave dominates for Ep < 10" 8 a.u.. 

The total annihilation cross-section is shown on Fig. || It follows the 1/v-law for Ep < 10~ 8 a.u. and decreases 
nonmonotonously for Ep > 10~ 8 a.u. Such nonmonotonic behavior is originated by the contribution of nonzero angular 
momentum partial waves, which is explicitly seen on Fig. m It is interesting to compare this cross-section with a 



semi-classical calculation 11 obtained under the following assumptions: i) the p motion can be treated classically 
and ii) the annihilation takes place with unit probability as soon as the p — H distance is smaller than a critical 
radius R c — 0.64re. The semi-classical cross-section, shown on Fig. is approximately 2.5 times bigger than our 
values for Ep < 10 -8 a.u. . This indicates that the low-energy pH annihilation is sensitive to the quantum dynamics 
of Protonium formation and could hardly be reproduced with models in which the details of such dynamics are not 
taken into account. 

The population of different Protonium states, calculated for energies Ep < 10~ 8 a.u. is shown on Fig. Protonium 
is produced primary in the S-states with principal quantum number 26 < n < 30. The P-states population does not 
exceed 15% of the whole captured fraction of antiprotons. These results confirm our qualitative estimation concerning 
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the channels which give the main contribution into the reaction amplitude in the low-energy limit. It is worth to 
mention that Protonium S-states population dominates only for p energies less than 10~ 8 a.u., while the population 
of states with higher I should increase with increasing energy |L7| . 

We conclude this paragraph by emphasizing that the Hp scattering observables significantly change their behavior 
at Ep ~ 1CP 8 a.u., a characteristic energy for the reaction. We will demonstrate that this behavior is determined by 
the presence of nearthreshold Hp bound and virtual states generated by the polarization potential. 

B. Nearthreshold metastable states 

The polarization potential is known to significantly modify the low energy cross-sections of atomic reactions. It 
plays an essential role in the Hp scattering. This potential produces a rich spectrum of Hp weakly-bound and virtual 
states which results from the long range character of the polarization forces and the heavy (in atomic scale) 
antiproton mass. Such states, being nearthreshold S-matrix singularities, determine the energy dependence of the Hp 
scattering cross-section. The main properties of such states and their relation with the observables are discussed in 
this subsection. 

We first remark that the polarization potential V po i alone generates several p weakly-bound states. The energy 
levels and mean radii of several nearest to the threshold S-states produced by V po i alone are shown in Table [j] (values 
marked by subscript II). These states are extremely prolonged and have very small binding energies. By switching on 
the short range part of the interaction, i.e. the complex nonlocal effective potential V e ff and the screened Coulomb 
V cs , the spectrum is modified and inelastic widths appear. Nevertheless, the main features, small binding energy 
(1CT 8 < Ebound < 10~ 3 a.u.) and big radius (4 < x < 27) r_g, remain. 

In the threshold vicinity the elastic S-matrix for L=0 is dominated by its singularities and can be written in the 
form: 

™*>-nj£t} w 

i 1 1 

where Z{ are the S-matrix poles with Re(zi) < due to Im(V e ff) < 0. In Fig. |?] are shown the trajectories of several 
S-matrix poles (Pi) and corresponding zeros (Zi) as a function of the strength of the V e ff imaginary part. As it can 
be seen, the presence of the negative imaginary part in the effective potential results in shifting the S-matrix zeros to 
the right into the IV and I quadrants, with the corresponding symmetrical shift of S-matrix poles into the II and III 
ones. The position of the nearest to the origin S-matrix zero (and pole) corresponds to an energy of E c ~ 10~ 8 a.u. 
and plays a role of characteristic energy for the reaction ([!]) . We notice however that this nearest to the threshold 
S-matrix singularity lies on the non-physical sheet, i.e. Re(k)< and Im(k)< 0, and corresponds to a virtual state. 
Its wave function has an exponentially incresing asymptotic and does not represent a physical state. 

As far as the usual definition of effective range can not be applied to the 1/R 4 polarization potential jL4|, we 
introduce the characteristic range of Hp interaction as: Ra — l/|feo|j where fco corresponds to the position of the 
nearest to the threshold S-matrix singularity. One can see, that for k > ko the scattering length approximation is no 
longer valid and higher order terms in the scattering amplitude expansion should be taken into account. With the 
result on Table I one gets Ra ~ 103 a.u. . 

It is seen from (^6|) that for antiproton incident energies Ep -C I0 -8 a.u., \k + Zi\ \k — Zi\ and so |5| — * 1. This 
explains why the inelasticity S% for Ep <C 10~ 8 a.u. turns to be much less than unit. For Ep > 10~ 8 a.u., and because 
there are several S-matrix zeros situated to the right from —z\, one has |fc + Zi\ < \k — and the reaction probability 
increases. 

To illustrate how the position of the nearthreshold S-matrix singularities determines the low energy scattering, 
we have calculated the inelasticity as a function of the dipole polarizability ay for a fixed energy (Ep — 10~ 8 and 
E = 10~ 6 a.u.). This function is plotted on Fig. ^. The strong oscillations between its maximum and minimum 
values with decreasing ad correspond to the motion of an S-matrix pole from the II to the III k-plane quadrant, while 
the symmetric S-matrix zero moves from the IV to the I quadrant. This means, that a weakly bound state becomes 
virtual. As it is seen from (^6|), the inelasticity reaches its maximum value when an S-matrix zero crosses the real 
k-axis. 

This last result shows that sufficiently high accuracy of calculations is required to obtain the scattering length 
value. In the same time the reaction amplitude for energies Ep 3> 10~ 8 a.u. is less sensitive to the exact position of 
the nearthreshold singularities and can be more easily calculated. We estimate our accuracy in the scattering length 
results to be about 30%. This uncertainty appears mainly from approximation used for V po i at short distances. To get 
more precise results one should increase the number of accounted closed channels and take into account the difference 
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between p and R in (|3[). Such corrections seem not to be important for understanding the physics of the treated 
process and are beyond the scope of present paper. 

We would like to emphasize that the nearthreshold character of mentioned S-matrix poles and zeros is determined 
by the long-range polarization potential. In the same time their exact position in complex k-plane can not be 
obtained without a proper treatment of the Protonium formation dynamics. In particular, the semiclassical black 
sphere condition does not hold in the energy domain Ep < 10~ 6 a.u. . In terms of S-matrix analytical properties the 
coupling with Protonium production channels produces comparatively big (for the energy domain of interest) shifts 
of the real part of the S-matrix zeros and reduces the inelasticity. 



C. Local approximation of the effective potential 



It was shown that the energy dependence of the reaction probability is determined by the existence of several 
nearthreshold states generated mainly by the long range polarization forces. This suggests the possibility to obtain 
a local complex potential which would be equivalent to the full Hp interaction in the energy range of interest. By 
equivalent we mean not only to reproduce the same reaction probabilities but to support the same nearthreshold 
spectral structure as well. 

We search for such equivalent local complex potential as a sum of three different terms 

Vi oc (R) = V S (R) + V CS (R) + V pol (R) 

V CS (R) and V po i(R) being respectively the Coulomb screened and polarization potential used in the previous section 
and V s a local short rang part to be determined. It was assumed to have the form: 

rr, m j -V ie - ai (^) -lW ie ~ Pl ^) R<R C ( ^ 

Vs{R) = < r ( r\ 

{ -iW 2 e P2 UbJ if R>R C 

and a satisfactory fit is obtained with the following parameter values: V\ = 0.572, W\ — W2 = 0.040, ol\ = 1.20, 
/3 1 =0 2 = 3.20 and R c = 2r B . 

In Table || the results of calculations in nonlocal effective potential and mentioned above local approximation are 
compared. They agree within few percent accuracy in the energy range 0.5 10 -9 — 0.5 10~ 6 a.u. 



D. Hydrogen-Antihydrogen interaction 

The results obtained for Hp interaction can be used for a qualitative treatment of different atom-antiproton (Ap) 
and atom-antiatom (AA) system. 

It is of particular interest to estimate the HH annihilation cross-section, and thus to examine the reaction: 

H + H -> Pn* + (e+e - ) (28) 

The HH system interacts at long distances via a dipole-dipole potential Vdd ~ — 6.5/i? 6 . This potential also generates 
a spectrum of nearthreshold states. Some of them with L=0 are shown in Table III. In analogy with the (Q) case one 



can expect that the corresponding S-matrix singularities will determine the (|28|) reaction dynamics. 

A qualitative estimation of the HH potential can be obtained by adding to the same short-range part as in Hp case 
the dipole-dipole long-range tail Vdd- The reaction (p8j) cross-section calculated in such a way is shown on Fig. 0. 
The characteristic energy for this reaction was found to be ~ 10~ J a.u. , corresponding to the position of the nearest 
to the threshold S-matrix singularity (virtual state with energy —7.8 10~ 6 a.u.) 

A similar treatment can be used to estimate the inelasticity energy dependence for different Ap or AA systems in 
the low energy limit. For such a purpose one has to find the nearest to the threshold S-matrix singularity, generated by 
polarization potential. The necessary condition for the validity of such a qualitative approach is that the characteristic 
range Ra of the Ap or AA long-range interaction should be much greater than the inelastic range ta- 

As it is seen from ( p4] ) the inelastic range is mainly determined by the mean radius of the last Protonium state 
open channel and thus given by: 

M A _ 

2n 2 1 

Ta ~ = — 

M A I a 
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where M A and I a are the Ap reduced mass and the first ionization potential respectively, n is the principal quantum 
number of the last open channel. A similar estimation for the AA inelastic range r AA can be obtained, if we take into 
account that Positronium is produced in this collision: 



M 



AA 



2n 2 



21 a - ep s 

Ih/2 
2n 2 



rA 



M AA 



1 



(21 A - I H /2) 



Here M a a is the reduced mass of the AA system, Ep s is the Positronium ground state energy. 

Like in the Hp case, the presence of nearthreshold virtual states may considerably increase the characteristic range 
of Ap or AA interaction. However it can be interesting to have a simple approximation of this range in the aim 
of comparison with preliminary estimations (see also fll2[). This is provided by the semiclassical condition for the 
number N of states : 



2M A V p A ol (R)dR 



ttN 



This condition may be rewritten as follows: 



For the Ap case we obtain: 



Ra 



nN 



Ra 



V2M A Cj 
^2M A C£/L(L 



if L = 
if L > 



(29) 



while for A A one has: 



Ra 



V2M aA C£ if L = 
i/2M AA cf/L(L + 1) if L > 



(30) 



C A and C A being the atom charge-dipole and dipole-dipole van der Waals constants. Finally, we get the following 
ratio of inelastic and polarization range: 



TA_ 

Ra 



I A ^ JlM A C{ 



I A y/lM A C{ 



if L = 



if L > 



(31) 



for atom-antiproton, and: 




i 77 T if L = 

(2/Wff /2) $/2M AA Cg 



(32) 



for atom-antiatom interaction. 

The ratios (|3l||3^), calculated for a wide range of different atoms, turns to be much smaller than unit in case L=0. 
In particular, for He, the less polarizable atom, they are ~ 0.02 for Hep, and ~ 0.05 for HeHe. The polarization range 
dominates over the inelastic one in the partial waves up to L ~ 10 for Hep, and L ~ 4 for HeHe. These values of L 
characterize the maximum angular momentum, which makes possible the existence of extended polarization states. 



IV. CONCLUSION 



A coupled channels model describing the Hp system at energies less than 10 6 a.u. has been developed. The results 
thus obtained substantially differ from the low energy extrapolations of the black sphere model and other classical or 
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semiclassical approaches. They show that such a low energy requires a quantum mechanical treatment in which the 
dynamics of the Protonium formation is properly taken into account. 

The effective Hp optical potential has been calculated in the framework of the coupled channels model. In this 
framework, the Hp scattering length and zero energy elastic cross section were found to be a — (—7.8 — ill.5)rg and 
a e i = 2426. 4r^ respectively. The Hp inelastic cross section has been calculated in the energy range from 10~ 9 to 
1CP 6 a.u.. It follows the 1/v behavior up to energies ~ 10 -8 a.u. where the scattering length approximation is valid. 
The inelasticity turned to be much smaller than the black sphere model predictions. 

The Protonium formation spectrum for the energies less than 10~ 8 a.u. has been calculated. We have shown 
that the population of S-states with principal quantum number from 26 to 30 acounts for 75% of the total captured 
fraction. 

The reaction dynamics is found to be determined by the existence of several nearthreshold states. Such states are 
produced by the long-range polarization potential and are shifted in the complex momentum plane by the coupling 
with Protonium formation channels. The Hp scattering length appears to be very sensitive to the position of the 
mentioned singularities and requires accurate calculations. 

A local approximation of the effective potential has been proposed for further applications. It reproduces the 
scattering observables in the considered energy range and has the same nearthreshold spectral properties. 

A qualitative extension of this approach to more general systems (atom-p and atom-antiatom) has been discussed. 

The results discussed in this work have been obtained within an approximate model. In view of them and motivated 
by the futur project of storing antimater at CERN it would be interesting to check the validity of the different 
approximations by developing more accurate treatments including an exact solution of the three body problem. 

V. ACKNOWLEDGMENTS 

The authors would like to thank I.S. Shapiro for suggesting the problem. One of the authors (A.V.) would like to 
thank D. Morgan for useful discussions. 

APPENDIX A: APPENDIX 

The aim of this appendix is to find the dominant channels in the expansion ( JTo| ) of the three-body wavefunction. 
We first analyze the behavior of the component $2 a t distances R 3> Tb- The equations system for $2 an d m 
terms of the projection operators P (0) and F reads: (|l9|): 

(H pp + PW ep P - Ep)^) + PW ep (l - P)F\^ d 2 ) + 

+ PW e p(l-P)(l-F)m)=0 (Al) 

(H ep + H pp + P(l - P)W ep (l - P)F - E) |$2*) + 
+P(1 - P)W ep {\ - P)(l - F)m + (1 - P)FW ep P\<P 1 ) = (A2) 



H ep + H pp + (1 - F)(l - P)W ep {\ - P)(l — F) — EJ m) + 

+ (1 - F)(l - P)W ep PF\^ + (1 - P)(l - P)VK e pP|$!> = (A3) 



By taking into account that at big R the projection operator F vanishes, equations (A1-A3) simplify into the 
following system, valid for R^> rs- 

(H pp + PW ep P - £p) |$i) + PW ef (l - P)m) - (A4) 
(Hpp + H ep + (1 - P)T^ e p(l -P)-E) \<P p 2 ) + (1 - P)VK eF P|$i) = (A5) 
This system can be solved with respect to $2 an d gives for the component $1: 

(h pp + PW ep P - Pp) |*0 + PW ep {\ - P)G po; (l - P)I^ e pP|$i> - (A6) 
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in which 



G po i = -(H pP + H ep + (l-P)W e¥ {l-P)-E) 1 (A7) 



The last term in equation (A6) is the polarization long-range interaction: 

V P oi = (<j>u\W e p(l - P)G pol {\ - P)W ep \4> ls ) (A8) 
The asymptotics of the Green-function G po i at R,R' 3> tb is: 

Gpo^R'.r.r 1 ) = £ J- (e"^ - S Q e~^ +ii ')) Q (r)0 Q (r') 



where p a = y // 2M(\E — e a \) and a is a set of spherical Coulomb quantum numbers. If we take into account that at 
big R, Vpoi acts on the very slowly changin g fu nction x{P) (the oscillation period of \ f° r P ^ r B is indeed much 
greater than rg) we can substitute G po i in (|Aq) by the following expression: 

G vol {R,R',ry) = 5{R-R') ^ { : )Mr,) 

Sli s a 

a 

By keeping terms up to 1/i? 4 we obtain the well-known charge-dipole potential asymptotic behavior: 

V pol (R^r B ) = -^ 



ad 



2 V (cp ls \d\<f>a) {<jy a Wis) (A9) 



sb - e a 



Here d stands for the dipole momentum operator. One recognizes in ( |A9| ) the expression for the Hydrogen dipole 
polarizability, ad = 9/2. We can thus conclude that the contribution of $2 at big distances R can be taken into 
account by introducing in the clastic channel the polarization charge-dipole potential. 

To qualitatively estimate the contribution of different channels at distances R ~ tb we first obtain the solution of 



(Al) in the distorted wave approximation. Component $5 S , obtained by neglecting the coupling to other components, 



is &i = 4>isX° with x° satisfying the equation: 
The contributions of components $2 and $2 are characterized by the integrals: 



(-^d 2 R + V cs (R)-E + e B )x° = 



(x°P\W ep (l~P)\F 

(A10) 

(x°P\W ep (l-P)\(l-F) 



An estimation of integrals (AlC) can be obtained if we take into account the semiclassical character of the wave- 
function x° and the Coulomb wavefunctions at R w tb in expansions (|l7| ) and (|l8|). We are dealing with an integral of 
fast-oscillating functions which has significant values only if there exist stationary phase points inside the integration 
region. The equation for such stationary phase points is: 

R J R 2n 2 

It can be shown that there are no stationary phase points for $31 while the contribution of $9_at the distance Rk, rs 
is mainly exhausted by Protonium states with principal quantum number 26 < n < 40 (see |l5| , |l6| ). 

We have, in conclusion, that in the energy domain of interest and large internucleon distances the c ontr ibution of 
$2 is the only important one and can be taken into account by introducing the polarization potential ( |A8| ) , while at 
the distances R ~ re the component $2 dominates, and can be described by a limited number of channels. These 
qualitative arguments are important for construction of the first approximation, and should be proved by further 
numerical calculations. 
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-2.5 10 


-6 _ 


i .2 10~ 7 


27.0 


-3.6 10 


-5 


-7.0 10- 


5 


i 8.4 10~ 6 


11.3 


-2.6 10 


-4 


-4.1 10- 


4 _ 


i 3.2 10~ 5 


7.3 


-9.2 10 


-4 


-1.5 10" 


3 _ 


i 8.6 10~ 5 


5.3 


-2.3 10 


-3 


-4.2 10" 


3 _ 


i 2.0 10- 4 


4.2 



TABLE I. Energies, Auger widths and mean radii (a.u) of L=0 Hp states. We denote by index I the results in the V po i alone 
and by index II those obtained with the full interaction (V po i + V cs + Veff) 



E p (a.u.) 


S A r (I) 


S'i (II) 


S (I) 


s (II) 




0.5 10" 9 


0.043 


0.043 


0.978 + i0.01 


0.978 + ^0.014 




0.5 10~ 8 


0.12 


0.122 


0.937 + xO.013 


0.936 + i0.021 




0.5 10~ 7 


0.266 


0.266 


0.836 — zO.185 


0.836-^0.177 




0.5 10- 6 


0.42 


0.425 


0.023 - i0.756 


0.034-^0.757 





TABLE II. Inelasticity (Sr) and S-matrix (S) values calculated in the full effective potential (index I) and in its local 
approximation (index II) at different energies (Ep) 



Ei 


E n 


Xn 


-7.8 10 


-b 


-6.1 10 


- b + i 


1.8 10- 


-b 




-1.9 10 


-4 


-4.3 10 




2.2 10- 


-4 


4.6 


-2.9 10 


-3 


-5.2 10 




1.2 10- 


-3 


2.8 


-1.1 10 


-2 


-2.9 10 


- 2 -i 


8.4 10- 


-3 


1.5 


-3.3 10 


-2 


-5.8 10 




9.2 10- 


-3 


1.3 



TABLE III. Energies, Auger widths and mean radii (a.u) of L=0 HH states. We denote by index I the results in the V po i 
alone and by index II those obtained with the full interaction (V po i + V cs + Veff) 
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FIG. 1. Imaginary part of effective potential V e ff(R,R' = R) 
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FIG. 2. Imaginary part of effective potential V e ff(R = 0.5, R') 
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FIG. 3. Real part of effective potential V e ff(R = 0.5, R') 
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FIG. 4. Inelasticity S r = 1 - |S| 2 for the H +p -> Pn + e reaction 
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FIG. 5. Annihilation cross section for H + p — * Pn* + e 
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FIG. 6. Protonium formation probabilities in states with different quantum numbers 
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FIG. 7. S-matrix nearthreshold zeros (Zi) and poles (Pi) 
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FIG. 8. Inelasticity for reaction H + p — > Pn* + e as a function of the dipole polarizability ad 
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